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We calculate the dominant contributions to the muon g — 2 at the two-loop level 
due to a generic pseudoscalar boson that may exist in any exotic Higgs sector in 
most extensions of the standard model. The leading effect comes from diagrams 
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f— ^ ' of the Barr-Zee type. A sufficiently light pseudoscalar Higgs boson can give rise 

f— ^ . to contribution as large as the electroweak contribution which is measurable in the 

1-^ ' next round of g — 2 experiment. Through the contribution we calculated here, the 

-L . anticipated improved data in the recent future on the muon g — 2 can put the best 

(— I I limit on the possible existence of a light pseudoscalar boson in physics beyond the 



standard model. 



PACS numbers: 14.80. Cp, 14.60.Ef, 13.40.Em 



^Work supported by the Department of Energy, Contract DE-AC03-76SF00515 



Precision measurement of the anomalous magnetic moment of the muon, a^ = ^{g^ — 
2), can provide not only a sensitive test of quantum loop effects in the electroweak 
standard model (SM), but can also probe the effects of some potential "new physics". 
The experimental average in 1998 Particle Data Book gives[T|] a*fP = 11659230(84) x 
10^^'^(±7.2 ppm). Recent measurements by E821 experiment at Brookhaven give 
a;fP = 11659250(150) x 10-i°(±13 ppm) (1997 data) and a^^P = 11659191(59) x 10"^° 
(±5 ppm) (1998 data). 

The E821 experiment has just announced the most accurate result]^ from the 1999 
data sample. The weighted mean measurement is 

a;fP = 116592023(151) x 10~^^(±1.3 ppm) . (1) 

Its high precision poses a direct challenge to the theoretical prediction of a^, which are 
usually divided into four sources, 

«/. = af ° + a?'''""" + «r + ^«M, (2) 

representing QED, hadronic, electroweak and the exotic (beyond the standard model) 
contributions respectively. The QED loop contributions have been computed to very 
high order [Q 

OED = ^^ 0.765857376(27) (-) + 24.05050898(44) (-) 
27r \7iJ Vvr/ 

+ 126.07(41) ('-') + 930(170) T-] . (3) 

The most precise value for the fine structure constant a = 1/137.03599958(52) can be 
obtained by inverting the similar formula for the electron Qe — '^ from the data[^. This 
gives 

^QED ^ 116584705.7(2.9) x 10"^^ , (4) 

with precision much higher than the expected experimental reach. The SM electroweak 
contribution up to two-loop level gives af^ = 152(4) x 10~^^||, for sin^ Ow = 0.224 and 
Mh = 250 GeV (In comparison, the one-loop SM electroweak contribution is 195 x 10~^^). 
The hadronic contribution]^ due to the hadronic vacuum polarization diagram has the 
largest uncertainty from the strong interaction, aHadromc _ 5739(57) x 10^^^. However, 



such uncertainty is still smaller than the experimental error in a^, even before results from 
future planned experiments which intend to measure the hadronic vacuum polarization 
directly. The total value in standard model is|0, 

a^(SM) = aQ^° + aj^^'^''°'^^= + af^ = 116591597(67) x 10-i^(±0.57 ppm). (5) 

The current measurement is well above the standard model prediction by a 2. 60" effect, 
Aa^ = (42.6 ± 16.5) x 10"^°. Given that a^'"^'"''''' and af^ are both positive, one can 
conclude that the current data have already probed these contributions, but they are not 
enough to fit the data. Extra contributions from new physics is needed. At 90%CL, Aa^ 

liesP] between (+21.5 1-63.7) x 10^"*^°. Note that the negative a|^^ (two-loop) plays an 

important role in claiming this discrepancy. 

Even without the recent experimental improvement, g — 2 data has already provided 
non-trivial constraints [T^ on physics beyond the standard model. For example, the con- 
straint on the minimal supersymmetric standard model (MSSM) due to its one-loop con- 
tribution to g — 2, via smuon-smuon-neutralino and chargino-chargino-sneutrino loops, 
is well known |TT[]. The resulting constraint depends on the masses of supersymmetric 



particles and tan /3. 

In theories beyond the standard model, there are usually additional scalar or pseu- 
doscalar bosons. In particular, some of the pseudoscalar bosons can potentially be light 
because of its pseudo-Goldstone nature, accidentally or otherwise. However, in collider 
search, it is known that searching for the pseudoscalar neutral boson is much harder than 
the scalar neutral or charged one. Therefore it is particularly interesting to see if one can 
constrain the pseudoscalar boson using low energy precision experiments. In this paper 
we wish to report that if the extended theory has a light enough pseudoscalar boson, its 
two- loop contribution to muon g — 2 can be as large as the one- loop electroweak effect. 
As a result the muon g — 2 can provide a very strong probe on a large class of theories 
beyond the standard model. 

The one-loop contributions to g — 2 from a scalar or pseudoscalar boson have been 
presented many times in the literature [|1^]. Besides two powers of m^ that is required by 
kinematics and definitions, the one-loop contribution is further suppressed by another two 
powers oirrin/Ma. However, the result is enhanced by a logarithmic loop factor, Inm^/Ma, 
coming from the diagram in which the photon is emitted by the internal muon. Therefore 



for a light enough Higgs mass, some hmit can be derived from g — 2 data just based on 
one-loop result. Nevertheless, as we shall see later, the two loop contribution is typically 
larger than the one-loop one by a factor of 2-10 for Higgs mass from 10 - 100 GeV. In 
addition, the one-loop and two-loop contributions have different signs for both scalar and 
pseudoscalar contributions. Therefore the one-loop contribution actually partially cancels 
the larger the two-loop contribution. 

The two-loop contribution of a scalar boson has been calculated in Ref. P, |^ in the 
context of the standard model. The contribution of any scalar boson beyond the stan- 
dard model can in principle be extracted from that calculation and we shall not dwell 
on this here except to note that the scalar boson gives negative contribution while the 
pseudoscalar gives positive contribution to Aa^. Also, we have parameterized our input 
Lagrangian as model independent as possible in order to make our gauge invariant result 
widely applicable to a large class of models. 

For Higgs mass larger than roughly 3 GeV, the dominant Higgs related contribution to 
the muon anomalous magnetic moment is through the two-loop Barr-Zee type diagram||13||, 
as in Fig. 1. Compared with the one-loop graph, the Yukawa coupling of the heavy 
fermion / in the inner loop together with the mass insertion of the heavy fermion in 
the two-loop graph will give rise to {nif/m^Y enhancement which can overcome the 
extra loop suppression factor of a/lGvr^. The internal gauge boson can be a photon or a 
Z^ . The Z^ contribution is typically smaller by two order of magnitude. It is included 
in this manuscript just for completeness. Note that unless CP violation occurs in the 
Higgs potential, there is no two-loop Barr-Zee type contribution to g^ — 2 associated with 
pseudoscalar boson and an inner gauge boson loop. 
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Fig. 1 The dominant two- loop graph involving a pseudoscalar boson that contributes to gi — 2. 
The cross location denotes a possible mass insertion. 

The form of the gauge invariant vertex function T'^'^ of a pseudoscalar boson a° of mo- 
mentum (p) turning into two photons {—k, fi), {q, v) due to the internal fermion or gauge 
boson loop is 

r'^^ = P(g2)e'^^"Vg/3 • (6) 

In general, the heavy fermion generation dominates in the loop. The Yukawa coupling is 
parameterized in a model independent expression, 



^^ ' <" 87r2 Mw Jo m^-z{\ 



Integrating the fermion loop momentum, we obtain the form factor 

87r2Mvi/ k m) - z{\ - z)q^ 

where rnj and g/ are the mass and the charge of the internal fermion in the loop. The 
color trace gives A'^ = A*"* = 3, Nl = 1. The above vertex is further connected to the 
lepton propagator to produce anomalous magnetic dipole moment aj"" for the lepton i, 



Svr^sin^^H^ M^ /4/ ^^' 'm^ W) 



ri In ,f . dz 

^(^> = I T^W^) ■ (^°> 

jr(l) = ^Cl2(f), with the Clausen's function Ch{9) = - /q^ In (2 sin f ) d9. As x > 1, 
xT{x) has the asymptotic form 2 + Inx. On the other extreme limit x <^ 1, T{x) 
approaches to ^ + In^ x. Our result is consistent with that from an unphysical Higgs 
boson in SM[0]. 

For the graph with the inner photon replaced by Z^ boson, its contribution to a^ can 
be calculated in a similar fashion, 



zo^o ^ a^mlAf^gl. ^ NiAfqfg(.m} 

"^ 87r2 sin^ 9w cos^ 9wM^z f^b,r ^z - ^l 
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with Qy = ^T^ifi) — qfsm^9w Note that, for both pseudoscalar and scalar boson 
contributions, only the vector coupling of Z^ to heavy fermion contributes to the effective 
vertex due to Furry theorem. Numerically, this Z^ mediated contribution turns out to 



be about two order of magnitude smaller than that of the photon mediated one. One 
suppression factor comes from the massive Z^ propagator and the other one comes from 
the smallness of the leptonic vector coupling of Z° boson, which is proportional to (— | + 

sin^^VK) 0.02. 

Taking the pattern of Yukawa couplings in MSSM as an example, we set Aj as cot /5 
(tan/3) for the u (or (i)-type fermion. The contributions due to top quark t , bottom 
quark h and tau lepton r in the loop respectively as well as the total are displayed in 
Fig. 2 for both tan/? = 30 and 50. In this MSSM pattern the t contribution is insensitive 
to tan/?. In addition, both the h and the r contributions, which are roughly the same 
order of magnitude, dominate over that of the top quark one for large enough tan /? and 
light enough pseudoscalar mass M^. For M^ ~ 15 GeV, the r contribution is larger than 
the 6-quark contribution. The total two-loop photonic contribution from the pseudoscalar 
boson, d?^ , can be as large as 10^^ for a large tan/5 when Ma < 10 GeV as shown in 
Fig. 2. For example, for Ma = 10 GeV and tan/3 = 50, a^;'^" (2-loop) = 1.2 x 10~^ which 
is above the upper limit allowed by the current experiment bound. Generically, for Ma ~ 
80-100 GeV, tan/3 ~ 50, a^ ranges in (7 - 9)xl0^^*^ which is close to the electroweak 
contribution. Note that the pseudoscalar contribution has the same sign as the hadronic 
or electroweak contributions. To derive constraint from the data one must combine the 



one- and two- loop contributions. The well-known one-loop contribution [|12| due to the 
pseudoscalar a" is 

«?(l-loop) = -^?^ f^V^ f^l , with H{y) = r "'^\ . (12) 
^^ ^' M^Ma^\2Mw ) \Mlj ^^' Jo 1-x + x^y ^ ' 

For small y, H{y) — > — Iny — ^ > 0. Note that the one-loop contribution is always 

negative in contrast to the two-loop contribution. In Fig. 2, we draw the absolute value of 

the one-loop contribution for easy comparison. For small Ma such as 10 GeV, the one-loop 

contribution can be as large as half of the two loop contribution and produce cancelling 

effect in g^ — 2. Complete cancellation occurs around 3 GeV for large tan/?. However the 

one-loop effect becomes smaller for larger Ma due to its additional suppression factor of 

{mf^/M^) ln(M^/m^) and basically negligible for Ma > 50 GeV. 

The up-to-dated measurement from E821 (1999) has already indicated 2.6 a deviation 

from the standard model. The allowed contribution from new physics falls into a small 

interval, Aa^ between (-1-21.5 1-63.7) x 10^^° at 90%CL. The positive two-loop con- 
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tribution is able to fit the data for large tan/3 ~ 50 and Ma ~ 40 GeV, as illustrated in 
Fig. 2. Note that for Ma lighter than roughly 3 GeV, the negative one- loop contribution 
dominates and gives the overall negative Aa^, which is disfavored by the current E821 
data. 



In CP conserving MSSM, there is a lower bound [ll4|, piq] on Ma > 88 GeV, which is 
only based on partial scanning with certain choices of benchmarks in the MSSM. Further- 
more, in more general supersymmetric models or in general two or more Higgs doublet 
models^], very little can be said about the potentially light pseudoscalar Higgs boson. 



Experimental constraint [1^ on Ma from LEP data is correlated to a rather light scalar 
Higgs boson. The model independent nature of our calculation makes it possible to derive 
useful information of the pseudoscalar boson sector in any theory beyond the standard 
model using the hard earned data on muon g — 2. 

Note that in general multi-Higgs doublet models, the tan (3 factor in our analysis may 
be supplemented by additional factor of mixing matrix elements. In addition, in any 
specific model, there may be additional two-loop contributions, such as the ones involving 
the physical charged Higgs boson or the neutral scalar boson. We assume that these 
contributions do not accidentally cancel each other. Given that the experimental limit on 
the masses of the charged Higgs boson as well as the neutral scalar boson are already quite 
high, it is very unlikely they will cancel the contribution of a relatively light pseudoscalar 
boson. 

In conclusion, in this letter we report a set of analytic formulas for the two-loop 
contributions of a generic pseudoscalar boson to lepton anomalous magnetic moment. 
Such pseudoscalar bosons may exist in any theory beyond the standard model and they 
are typically harder to constrain using collider experimental data. In this paper, we show 
that strong constraint on such sector can be derived from the precision data on muon 
anomalous magnetic moment from the going and future experiments. We hope our work 
add importance and urgency to these low energy precision experiments. 
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Note added: While the original manuscript of this work was in the review process, the 
E821 experiment announced its new measured value. Based on this non-trivial result, we 
have updated our Fig. 2, which shows the implication on the mass of the pseudoscalar 
Higgs boson. 
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Fig. 2: The dotted lines plot the positive two- loop contributions from the inner t,b, t loops 
to g^ — 2 due to the pseudoscalar a^ versus Ma at tan/3 = 50, while the dashed-dotted line 
plots the negative contribution from the one-loop diagram. The sum in solid (dashed) curve 
shows cancellation at low Ma mass for tan/3 = 50 (30). The shaded areas shows the allowed 
contribution to Aa^ from new physics at 90%CL based on E821 measurement (1999). The 
region of negative Aa^, mainly from the one-loop contribution for a very light Ma < 2.9 GeV, 
is not favored by data. 
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